Triangular Trimers on the Triangular Lattice: an Exact Solution 
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A model is presented consisting of triangular trimers on the triangular lattice. In analogy to the 
dimer problem, these particles cover the lattice completely without overlap. The model has a 
honeycomb structure of hexagonal cells separated by rigid domain walls. The transfer matrix can 
be diagonalised by a Bethe Ansatz with two types of particles. This leads two an exact expression 
for the entropy on a two-dimensional subset of the parameter space. 

PACS number(s): 05.50.+q 

In the course of years a few exactly solvable lattice gas 
models have been found. The Ising model 0, proposed 
in 1920 by Lenz as a model of a ferromagnet, can 
be interpreted as a lattice gas with hard-core repulsion 
and short-range attraction. The (zero-field square lat- 
tice) Ising model was solved in 1944 by Onsager ||. It 
exhibits gas-liquid coexistence below a critical tempera- 
ture and a single fluid phase above. 

Another lattice gas is the hard hexagon model |Q, 
solved in 1980 by Baxter ||. It has a continuous fluid- 
solid transition. At high-density (solid) the particles se- 
lect one of three sub-lattices; at low-density (fluid) these 
sub-lattices are evenly occupied. 

As a final example we mention the dimer problem. It 
was solved for planar lattices in 1961, independently by 
Kasteleyn [|| and by Temperley and Fisher M . A dimer 
is a particle that occupies two adjacent lattice sites. As in 
the Ising and the hard hexagon model two particles can- 
not occupy the same lattice site. In contrast to these 
models it is also required that all sites are occupied. 
The configurations are coverings of a lattice with dimers, 
without empty sites or overlap. The dimer problem is 
reviewed in Ref . || . We discuss the dimer model on the 
honeycomb lattice in some detail, because it has illustra- 
tive similarities to a new model we shall introduce below. 

A configuration of the honeycomb lattice dimer model 
can be viewed as a number of domains consisting of ver- 
tical dimers, separated by zigzagging domain walls made 
up of dimers of the other two orientations. This is illus- 
trated in Fig. [lj. The domain walls run from the bot- 
tom to the top of the lattice, so that any horizontal line 
through the system meets all domain walls. Hence the 
number of domain walls is the same in each horizontal 
slice, in other words, it is a conserved quantity. 

Consider the entropy for fixed density p of domain wall 
dimers. From the exact solution of the model it can be 
calculated that for low p the entropy per dimer is given 
by 



FIG. 1. The dimer model on the honeycomb lattice. Do- 
mains of vertical dimers are separated by domain walls con- 
sisting of dimers of the other two orientations. To guide the 
eye the domain walls are shaded. 
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The linear term reflects the zigzag freedom of the do- 
main walls; each domain wall dimer contributes log 2 to 



the entropy. The cubic term is due to the (repulsive) in- 
teraction between the domain walls: when two domain 
walls meet some of the zigzag freedom is lost. 

Now give chemical potentials p to the dimers in the do- 
main walls and to the vertical dimers. For p < — log 2 
the free energy F = —pp — S(p) is is an increasing func- 
tion of p for small p, so no domain walls will be present. 
F° r A* ^ — 1°S 2 the free energy has a minimum at some 
small positive value of p. At p = — log 2 there is a transi- 
tion between a frozen phase consisting of vertical dimers 
only and a rough phase where dimers of all three orien- 
tations are present. 

Inspired by the dimer model we consider coverings of 
the triangular lattice by triangular trimers. A trimer is a 
particle that occupies three lattice sites. As in the dimer 
problem we require that there are no empty sites and that 
there is no overlap. Fig. ^| shows a typical configuration. 
In this paper we present our main results on this model. 
We intend to publish a more detailed account later. 

The model admits very regular configurations where 
the trimers occupy a sub-lattice of the triangular faces. 
There are six such sub-lattices, which we number 0, 1, 
. . . , 5 as indicated in Fig. ||. Note that the up and down 
triangles make up the even-numbered and odd-numbered 
sub-lattices, respectively. 

Consider configurations where trimers on sub-lattice 
predominate. They consist of hexagonal domains of 



1 





/ v° 


1 


2 


\ 3 / 


4 


\ 5 / 






4 







1 


2 


,3 /\ /\ /\ /\7 


/ \/° 


1 


2 


\ 3 / 




\ 5 / 


o\/ \/ \/ \/ \ 



FIG. 2. A typical configuration of the trimer model. Each 
lattice site belongs to precisely one trimer. We divide the 
lattice faces into six sub-lattices, numbered 0, 1, . . . , 5. Filling 
one sub-lattice completely and leaving the other five empty 
gives a very regular configuration of the model. 

trimers on this sub-lattice, separated by straight domain 
walls that form an irregular honeycomb network. There 
are three types of domain walls, of different orientations. 
The domain walls that run from lower right to upper left 
will be termed L; they are made up of trimers on sub- 
lattice 5. Those running from lower left to upper right 
will be called R; they consist of trimers on sub-lattice 1. 
The vertical domain walls are made up of trimers on sub- 
lattice 3. When domain walls of the three different types 
meet in a Y-shape a trimer on sub-lattice 2 or 4 occurs. 
This is illustrated in Fig. |^. 

For the time being we require that the model is 
isotropic, in the sense that there are equal amounts of 
the three types of domain walls. Let p denote the den- 
sity of domain wall trimers. The domain walls are rigid, 




FIG. 3. The trimer model has three types of domain walls, 
corresponding to the three odd-numbered sub-lattices. The 
trimers on sub-lattice are shaded dark grey; the trimers 
on the other sub-lattices are shaded light grey, and their 
sub-lattice number is given. 



so they have no zigzag freedom contributing to the en- 
tropy. Therefore the low density expansion of the entropy 
contains no term linear in p. There is, however, freedom 
in the sizes of the domains ||. For example, it is possible 
to enlarge a single domain while simultaneously shrinking 
its six neighbours. The contribution per domain depends 
on the linear dimensions of the domains, and is roughly 
proportional to — log p. The number of domains is ap- 
proximately proportional to p 2 . Hence the "breathing" 
entropy is given for low p by 

S « -Kp 2 log p, (1) 

where K is some (positive) proportionality constant. 

If a chemical potential p is given to the domain wall 
trimers, the free energy for low p is 

F w -fxp + Kp 2 log p. 

This is an increasing function of p for p < and a de- 
creasing function for p > 0. Hence the free energy takes 
its minimum either at p — or at a large value of p, 
for which the approximation (]]]) is not valid. For small 
p there are no domain walls, but when p passes some 
threshold p jumps to a positive value. Thus the phase 
transition is different from that in the honeycomb lattice 
dimer model, where the domain wall density increases 
gradually at the phase transition. 

Now we return to the model without the isotropy re- 
quirement. Let N denote the total number of trimers, Ni 
the number of trimers on sub-lattice i, and pi the partial 
density Ni/N. These six sub- lattice densities obviously 
satisfy 

PO +Pl + P2 +P3 + PA + P5 = 1- (2) 

It can be shown that, when toroidal boundary conditions 
are imposed, they also satisfy 

P0P2 + P2PA + PiPa = PiPi, + P3P5 + PbPi ■ (3) 

When the total density of down trimers p v = Pi + P3 + 
P5 is small, it follows easily from (||) and (|J) that one of 
Pa, P2 and p^, say po, is larger than the other two. If 
there is no further symmetry breaking po > Pi = P3 = 
P5 > P2 = Pi- By the same token when p v is close to 1 
the symmetry between the down sub-lattices is broken. 
Therefore when p v is increased from to 1 at least one 
phase transition is expected. 

Beside (§) and (§) we have found no more constraints 
on the sub-lattice densities. Therefore of the six sub- 
lattice densities four are independent. We would like 
to know the entropy (per trimer) as function of these 
four parameters. We have been able to compute it for a 
two-dimensional subset of the four-dimensional parame- 
ter space. The calculation is rather lengthy, so here we 
only give an outline of the method, and a description of 
the final result. 
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View each vertical domain wall as a combination of 
one L domain wall and one R domain wall. Then the L 
and R domain walls run without interruption from the 
bottom to the top of the lattice. Therefore the number 
of L domain walls and the number of R domain walls 
are constant throughout the system. This is analogous 
to the situation described above for the dimcr model on 
the honeycomb lattice, except that there are now two 
conserved quantities hl and Ur instead of a single one. 
We introduce the densities pl = nh/L and pr = njt/L, 
where 3L is the number of sites in a horizontal row of 
the lattices. They can be expressed in terms of the sub- 
lattice densities: 

Pl = 1 - Po - Pi + P3 + Pi, 

PR = 1 - PO + P2 + P3 - Pb- 

It is suggestive to interpret the vertical lattice direction 
as "time" and the horizontal direction as "space". The 
domain walls then are viewed as world lines of two types 
of particles, L and R. In a vertical domain wall one L and 
one R particle form a "bound state" . 

The model can be formulated in terms of a transfer 
matrix, which describes the "time" evolution of the sys- 
tem of L particles and R particles in one "space" dimen- 
sion. Solving the model boils down to determining the 
largest eigenvalue of this operator, or more precisely, its 
maximum over all particle numbers hl and ur. This 
we have achieved using coordinate Bethe Ansatz; the so- 
lution is similar to that of the square-triangle random 
tiling model, due to Widom [M and Kalugin pE. A 
model can be solved in this way only if, in some sense, 
the many-particle interactions factorise into two-particle 
interactions. It turns out that for the present model this 
is indeed the case. It is noteworthy that the L particles 
among each other are free fermions, as are the R par- 
ticles, but that the interaction between an L and an R 
particle is non-trivial. 

The Bethe Ansatz allows for numerical computations 
for the system on an infinitely long cylinder of finite 
circumference. These computations can be done to 
arbitrary precision, and effectively for the full four- 
dimensional parameter space of the model. In the ther- 
modynamic limit the Bethe Ansatz gives rise to a set of 
two coupled integral equations. The physical quantities 
we are interested in, such as the densities pl and pr and 
the entropy, can be expressed in terms of the functions 
satisfying these equations. These can be solved analyt- 
ically in a special case pT[] . Thus we have obtained an 
exact expression for the entropy for a two-dimensional 
family of sub-lattice densities. It will be seen below that 
this family is given by pi = p 3 = p 5 (or p = p 2 = p 4 ). 

This solution is parametrised by a complex number b 
with Im b > 0. Write 



(4) 





FIG. 4. The exact expression for the entropy is formulated 
in terms of contour integrals. The curves Cl and Dl in the 
right half plane are deformations of the line segment joining 
&l and 6 L , and Cl lies to the right of -Dl. No orientation of 
-Dl needs to be specified because it does not occur as an inte- 
gration contour but only as a branch cut. Analogous remarks 
apply to Cr and Dr. 



with Rc 6l > and Re 6r < 0. It follows that Im &l > 0, 
Im 6r > and 6l&r = — 1- Take contours Cl and Dl 
running from b L to 6l, and Cr and _Dr running from 6r 
to 6 R . The arrangement of these four curves must be as 
shown in Fig. ||, but their precise shape is immaterial. (In 
fact there are more solutions, corresponding to contour 
configurations different from that in Fig. ||. Since they 
are related by permutations of the six sub-lattices, we 
only treat one case here.) Define the complex function 



t(z) = 



- b 



1/6 



-b* 



Fix a branch ti,{z) with branch cuts Cr and Dl by 
*l(0) = exp(7ri/3), and a branch tji(z) with branch cuts 
Cl and £)r by tR(0) = exp(— 7ri/3). The domain wall 
density pl is given by 

Pl = — I dz. 
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and pr is given by the same equation with all subscripts 
L changed into R. The sub-lattice densities are 
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Define auxiliary integrals </>l and Sl by 

tL(z)+t L (zy 1 



Z -lb 



-b: 



dz, 
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The real parts of these integrals do not depend on the 
choice of the integration contours, which must not meet 
the branch cuts Cr and D^. The auxiliary integrals 0r 
and £r are defined analogously. The entropy per trimer 
is given by 



S = E L + Sr + -(2pR - + t (2pl - Pr)0r- 

6 o 

We started out with the model parametrised by the six 
sub-lattice densities satisfying the constraints (|^) and (||) . 
The exact solution described above has two parameters, 
so it covers a two-dimensional set of sub-lattice densities. 
The solution is however parametrised by a complex num- 
ber b, and not in terms of these densities. It follows from 
(||) that pi = ps = /?5 for the exact solution. The space 
of sub-lattice densities satisfying this constraint as well 
as (^) and ^ is two-dimensional. 

In the limit b — > the domain wall densities pl = PR 
vanish; the system is then filled with sub-lattice trimers 
and its entropy is zero. At b — 2i the system is in a sym- 
metric phase with all sub-lattice densities equal to 1/6; 
its entropy is 5 sym = log(3v3/4). When b is taken be- 
tween and 2i on the imaginary axis, &l and 6r lie on 
the unit circle. Then pl = Pr, so the sub-lattice densi- 
ties satisfy p2 = pi- For small p v < 1/2 this equation 
together with p\ — pz — ps describes the most symmet- 
ric case for the sub-lattice densities. Based on numeri- 
cal Bethe Ansatz calculations we believe that for given 
Pv < 1/2 the system takes its maximum entropy at these 
sub-lattice densities. Fig. || shows the entropy S as func- 
tion of p v - The entropy for p v > 1/2 was obtained from 
that for p v < 1/2 using the symmetry between the up 
and the down trimers. 

The entropy 5* is a convex function of p 7 for < 
Pv < 1/2- A system with p v m this interval is ther- 
modynamically unstable. It would separate into a phase 
with /9 V = and a phase with p v = 1/2, except for the 



fact that the model does not admit an interface between 
these two phases. Similarly a system with 1/2 < p v < 1 
would demix into phases with p v = 1/2 and p v = 1. 

The transition between these phases can also be con- 
trolled by assigning a chemical potential p to the down 
trimers instead of imposing their density p v - From 
Fig. U it is seen that for p < —2S sym the free energy 
F = — ppv — "5*(Pv) takes its minimum at p v = 0, 
so all trimers are on one of the up sub-lattices. For 
— 2iS sym < p < 25 sym the minimum of F is at p v = 1/2, 
so all sub-lattices are equally occupied. For p > 2S sym 
the minimum of F is at p v = 1 , so the system is again in 
a frozen phase. 
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